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Sine-Gordon field tlieory is used to investigate tlie phase diagram of a neutral Coulomb gas. A 
variational mean field free energy is constructed and the corresponding phase diagrams in two (2d) 
and three dimensions (3d) are obtained. When analyzed in terms of chemical potential, the Sine- 
Gordon theory predicts the phase diagram topologically identical with the Monte Carlo simulations 
and a recently developed Debye-Hiickel-Bjerrum (DHBj) theory. In 2d we find that the infinite order 
Kosterlitz-Thouless line terminates in a tricritical point, after which the metal-insulator transition 
becomes first order. However, when the transformation from chemical potential to the density is 
made the whole of the insulating phase is mapped onto zero density. 



PACS numbers: 64.70.-p, 05.70.Fh, 64.60.-i 

The Coulomb gas provides a paradigm for the study of 
various models of critical phenomena [y. In particular 
it is well known that the two dimensional (2d) Coulomb 
gas (CG) can be directly used to study the superfluidity 
transition in ^He films, arrays of Josephson junctions [^ , 
melting of two dimensional crystals [g| , roughening tran- 
sition Q, etc. Not withstanding its versatility our full 
understanding of the most basic model of Coulomb gas, 
namely an ensemble of hard spheres carrying either pos- 
itive or negative charges at their center, is still lacking. 

It is now well accepted that at low density the two 
dimensional plasma of equal number of positive and 
negative particles undergoes a Kosterlitz-Thouless (KT) 
metal insulator transition H. This transition is of an 
infinite order and is characterized by a diverging Debye 
screening length. Thus in the low temperature phase 
(insulator) all the positive and negative particles are as- 
sociated into the dipolar pairs, while in the high temper- 
ature phase (conductor) there exists a finite fraction of 
unassociated, free, charges. As the density of particles 
increases the validity of the KT theory becomes ques- 
tionable and the possibility of the KT transition being 
replaced by some kind of first order discontinuity has 
been speculated for a long time g. The idea that there 
can exist a discontinuous transition between the insulat- 
ing and conducting phases has gained further credence in 
view of the increasing computational power and an im- 
proving algorithm design needed for running large scale 
simulations of the particles interacting by a long range 
potentials ||7[]. Thus, it has been demonstrated quite con- 
vincingly that at high densities the KT infinite order line 
becomes unstable and is replaced by a first order coexis- 
tence between the low density insulating vapor and the 
high density conducting fiuid-like phase. From the theo- 
retical perspective, however, the nature of this metamor- 
phosis is far from clear. At the moment there are appear 
to exist two competing views of what happens to the two 
dimensional plasma at higher densities. The first of this 
theories, presented by Minnhagen (Mh) et al. B in an 



series of papers going back ten years, predicts that the 
KT line will terminate in a critical end point, while the 
critical point of the coexistence curve separating the low 
and the high density phases lies in the conducting region. 
It is important to note that within the Minnhagen's the- 
ory the portion of the coexistence curve in between the 
critical point and the critical end point has both vapor 
(low density phase) and the liquid (high density phase) 
which are conducting. Although the Minnhagen's ap- 
proach is often characterized as a version of Renormal- 
ization Group (RG), this is a misnomer. Its basis, which 
lies in a clever combination of a linear response formal- 
ism with some aspects of Sine-Gordon (SG) field theory, 
is much closer to the integral equations of liquid state 
theory than to the RG. The RG methodology being used 
more as a tool in studying the solutions of the integral 
equation found by Mh. 

An alternative approach suggested by Levin et al. H 
is based on a recently developed Debye-Hiickel-Bjerrum 
(DHBj) theory |jiq]. This method, which is intrinsically 
mean-field, relies on calculating the full electrostatic free 
energy of the ionic solution based on a linearized Poisson- 
Boltzmann equation. The effects of linearization are then 
corrected by allowing for the presence of dipolar pairs, 
the density of which is determined through the law of 
mass action. This theory has proven extremely powerful 
in elucidating the critical properties of a three dimen- 
sional (3d) electrolyte solutions iQ. In particular, the 
coexistence curve obtained on its basis was found to be in 
an excellent agreement with the recent MC simulations 
||ll| . The application of this theory to the 2d plasma 
has lead to a stark disagreement with the work of Mh. 
Where is Mh has found that the KT line terminates in 
a critical end point the DHBj theory predicts that it will 
terminate in a tricritical point after which the vapor in- 
sulating phase will coexist with a liquid conducting phase 
§ (see inset Fig.|). 

Since the DHBj theory is intrinsically mean-field, one 
might argue that the fluctuations, such as a variation 



in dipolar sizes, might modify the phase diagram. This, 
however, is not very hkely. It is weU known that a prop- 
erly constructed mean-field theory almost always retains 
the topology of the phase diagram upon inclusion of fluc- 
tuations. One of the few exceptions is when the volume 
of fluctuations is extremely large, such as in the case of 
transition between disordered and lamellar phases |I3] in 
magnets or amphiphilic systems. This, however, is not 
the case here. Furthermore the scaling of dipolar sizes 
can be included in a straight forward way into the DHBj 
theory, leaving the topology of the phase diagram un- 
changed |]l3|. The metal- insulator line then becomes in 
exact agreement with the KT theory, in particular, giving 
the correct critical exponent v = \/2 for the divergence of 
the screening length upon the approach to the transition. 
The tricritical point persists, while the first order coexis- 
tence curve remains extremely narrow in the vicinity of 
the tricritical point. 

Comparing the predictions of Mh and DHBj theory to 
MC simulations we find that neither one is in a quantita- 
tive agreement with MC, which finds that the first order 
transition appears at a temperature which is significantly 
lower, and the density which is significantly higher than 
the prediction of either one of the above theories . Nev- 
ertheless, the topology of the phase diagram observed on 
the basis of DHBj theory is the same as found in MC. Fur- 
thermore, the location of the tricritical point obtained in 
the MC corresponds closely to the region of the phase 
diagram where the narrow DHBj coexistence is found to 
swell significantly p^ ]. 

The current impasse lead us to reexamine some of the 
foundations on which our understanding of CG is based. 
Most of the rigorous theorems concerning the nature of 
interactions inside the neutral plasma are based on the 
isomorphism between the CG and the Sine-Gordon Field 
Theory |g|, (b)] ||lj|. The mapping is exact only for the 
point Coulomb gas in the grand canonical ensemble. The 
short range repulsion is included post facto by introduc- 
ing a suitable cutoff on all momentum space integrals. 
To what extend this procedure is valid is far from clear. 
Nevertheless, if the Sine-Gordon field theory is renormal- 
ized one obtains exactly the KT fiow equations in terms of 
renormalized fugacity and temperature |l5[ . This equa- 
tions, however, remain valid only for low density (small 
fugacity). 

The attempts to construct a Sine-Gordon based mean- 
field theory go back to the work of Saito ||l^, who has 
observed that already at the mean-field level the Sine- 
Gordon Hamiltonian knows about the metal insulator 
transition. In particular Saito was able to show that 
at low fugacity the Debye screening length diverged as 
^jj = e^/*", where t = {T-Tc)/Tc and v ^ 1. This should 
be compared with an equivalent expression obtained by 
KT but with I' = 1/2. Using field theoretic methodology 
Zhang et al. Il^ extended the mean-field type of calcula- 
tions of Saito and found that above a critical fugacity the 



screening length has a discontinuous jump from a finite to 
an infinite value. Zhang et al. then interpreted this point 
as a tricritical point terminating the continuous line of 
metal insulator transition. 

To compare the results of the Sine-Gordon based the- 
ory to the MC, one must be able to come up with a trans- 
formation from the fugacity, which is a natural variable 
in the field theoretical description, to the density, which 
is what the MC simulations measure. In the following 
we present a simple variational mean-field theory that 
accomplishes just that. It is in the process of transform- 
ing the phase diagram from the temperature-fugacity to 
temperature-density plane that the surprising new results 
were found. 

Our starting point is the Grand Canonical partition 
function for point particles of charge ±q 
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Here N = N^ + iV_ is the total number of particles im- 
mersed in a homogeneous medium of dielectric constant 
D and A = {h^ /27rmkBTy^^ is the thermal wavelength; 
the two-dimensional interaction term is U{rij) ~ Inry/a, 
where a is an arbitrary scale and Tjj is the distance be- 
tween particles i and j. The fugacity is related to the 
chemical potential through z± = e^^^ and, along with 
the temperature (/3 = l/fc^T), determines all thermody- 
namic characteristics of the two-dimensional Coulomb- 
gas (2D CG). 

To explore the thermodynamic properties of the above 
partition function it is convenient to map it onto the Sine- 
Gordon Field Theory [|, (b)] (l||. Thus the partition 
function Z can be expressed as a functional integral over 
a real field 6 
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is the effective Hamiltonian for a neutral Coulomb gas, 
and z = z{a/X)'^e^'^ u{o)/2D jg ^-^^ fugacity renormal- 
ized by a self-energy term. It is interesting to note that 
the saddle point of the Sine-Gordon Field Theory corre- 
sponds to the familiar Poisson-Boltzmann equation. In 



this paper, however, we will not use this analogy but 
instead construct a variational bound for the free en- 
ergy. To this end we shall rely on the Gibbs-Bogoliubov- 
Feynman inequality exploring the convexity of free en- 
ergy, Q < G = Go+ < H — Ho >o, where Go is the 
free energy associated with an arbitrary trial Hamilto- 
nian Hq. The brackets indicate averaging over Hq. It is 
particularly convenient to choose as a trial Haniiltonian 
one having a Gaussian form. 
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In this case the free energy Qo and the average < H 
Hq >o are easily calculated, and we find 
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To perform the momentum space integrals the ultraviolet 
cutoff (A = 1/a), corresponding to the effect of the hard- 
core, was introduced |18| . The optimum upper bound is 
found by minimizing the free energy G over all possible 



m . We find that the value of 



which leads 



to the optimal approximation to the real free energy G 
satisfies 
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where T* — ksTD /q^ is the reduced temperature. The 
parameter ttiq is inversely proportional to the Debye 
screening length, ^/j, inside the electrolyte solution, since 
it can easily be shown that the effective potential between 
two test particles separated by a distance r is Ve//(r) '^< 
(/)(O)0(r) >o. Furthermore if ttiq = (^_d = oo) signifies 
that there is no screening which means that all the ions 
have paired up forming dipolar pairs. If this is the case 
the presence of an insulating phase is assured. 

Indeed from the equations (|g) and (0) we find that 
the free energy possesses two minimas one of which is 
for TOg = and the other TOq ^ 0. The first order phase 
transition occurs when the free energies corresponding to 
the two local minimas become equal. The phase diagram 
in the {z — T) plane for the 2D Coulomb Gas is presented 
in Fig.|l|. It is essentially divided into two regions, each of 
which is characterized by a specific value of rn^: a con- 
ducting phase with a finite value of TOq and an insulating 
phase with TOq = 0. Separating these two phases there is 
a first-order transition line that ends at a tricritical point 
C {T* = 1/4 and z^ = l/167r). Below z^ and close to T*, 
the Eq. (0) can be approximated by 
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FIG. 1. Phase diagram of the 2D CG in the z — T* plane. 
The solid line corresponds to the first-order transition, while 
the dashed line is the KT infinite order metal-insulator tran- 
sition. The tricritical point C is located at z — l/lGvr and 
T* = 1/4. 

where t = I- 1/AT* and i^ = 1. When T* -^ 1/4 the 
equation defines a line of a critical points {^d — oo) that 
separates the conducting and insulating phases. This cor- 
responds to the usual KT line of metal-insulator transi- 
tions. Contrary to the appearance, the first-order line 
and the critical line join smoothly at the tricritical point, 
with the tangency of the first-order line ensured by the 
divergence of dz/dT* ~ In |t|/47r, when t ^ 0". 

As was emphasized in the introduction, in order to 
compare the results of our variational treatment with 
those of MC simulations it is essential to perform a trans- 
formation from the fugacity-temperature plane into the 
density-temperature domain. To this end we note that 
G is related to pressure and volume through G — —PV, 
while the density is p = zd{(3P)/dz. The transformation 
is then easily achieved and we find the coexistence curve 
presented in Fig||. This curve is topologically identical 
to that obtained on the basis of pure linearized Debye- 
Hiickel (DH) theory ||^. In particular we find that the 
high density conducting phase coexists with the zero den- 
sity insulating phase. Namely, although the SG theory 
knows about the metal-insulator transition it can not give 
a proper account of the low density phase. Instead of pro- 
ducing dipoles the oppositely charged ions self annihilate 
on contact! 

From our treatment it is not clear if this is a true prop- 
erty of the SG model or is an artifact of mean-field treat- 
ment, or maybe a result of the artificial way in which the 
hardcore was introduced into the model. Whichever the 
case, it is interesting to compare this result with the re- 
cently developed Debye-Hiickel-Bjerrum theory (DHBj) 
which, although also is of a mean-field type, does pre- 
dict a finite density for the insulating phase. The coexis- 
tence curve for the DHBj theory is presented in the inset 
ofFig|. 

The inability of SG model (at least at the mean-field) 
level to give a proper account of dipoles is also confirmed 
in rf = 3. In this case a first-order phase transition at low 
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FIG. 2. Coexistence curves for the 2d SG field theory. The 
tricritical point is at T^ =1/4 and p* = 0. The inset repre- 
sents the phase diagram for d = 2, according to DHBj theory 
[9] , with the KT line ending at a tricritical point (dashed line) 
locahzed at p* ~ 0.003954 and T^ = 1/4. 

temperatures between a low density (vapor) and a high 
density (liquid) phases is found as expected. Indeed, 
based on the DHBj theoretical [|l^ studies and simula- 
tions ||l^, the critical point is localized at T* ~ 0.057 
and p* ~ 0.025. The resulting coexistence curve that 
emerges from our variational treatment predicts a criti- 
cal point at T* ~ 0.0565 and p* ^ 0.00135. While the 
critical temperature T* is in agreement with previous re- 
sults, the critical density p* is too small. This values 
should once again be compared with the pure linearized 
DH theory which does not account for the existence of 
dipoles; in that case it was found that T* = 0.0625 and 
p* = 1/64tt ~ 0.005. The underestimate of the critical 
density clearly indicates that at least at the mean-field 
level the SG theory, just as pure DH theory, does not 
give a proper account of non-linear effects such as the 
formation of dipoles. 

How can these non-linearities be included is far from 
obvious. Why should the SG theory predict a metal- in- 
sulator transition in the temperature-fugacity plane only 
to later map the whole of the insulating phase onto zero 
density? What is the proper class of diagrams that would 
have to be summed to produce a finite density for the in- 
sulating phase? This questions require serious attention 
if we wish to have a complete theory. Inability of the 
Sine-Gordon model, at least at the mean field level, to 
give a proper account of the low density phase might 
also be responsible for the distinct predictions between 
Mh and DHBj theories. In particular, it can be shown 
that the variational method that we have used corre- 
sponds to the summation to all orders of a certain class 
of diagrams. In the case of the standard scalar field the- 
ory this is the familiar Hartree-Fock approximation [^ . 
This class of diagrams is obviously insufficient if we are 
to believe that the SG theory can give a realistic ac- 
count for the phase structure of the Coulomb gas. In 
his approach Minnhagen also relied on the SG theory 
to calculate the charge-charge correlation function. To 



this end he summed another set of diagrams. If that set 
was incomplete it could lead to some undesirable effects 
such as, for example, the wrong topology of the phase 
diagram. At the moment, however, this is only a specu- 
lation and a renewed theoretical effort is needed to study 
the Sine-Gordon field theory, now that it is evident that 
this model, besides the KT transition, also contains a 
first order discontinuity. 
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